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Abstract. We prove that linear instability implies non-linear instabil- 
ity in the energy norm for the critically dissipative quasi-geostrophic 
equation. 



1. Introduction 

A fundamental equation in oceanography and meteorology is the 3 dimen- 
sional Navier-Stokes equation in the context of a rapidly rotating, density 
stratified, viscous, incompressible fluid. Both the forces of rotation and 
stratification impose a tendency toward 2 dimensionality on the 3 dimen- 
sional fluid motion, and this leads to approximate and simpler mathematical 
models. Important non-dimensional parameters are the Ekman number (the 
strength of the viscous term relative to rotation) and the Rossby number 
(the strength of the nonlinearity relative to rotation). In many geophysical 
problems these parameters are very small. A set of approximations based 
on asymptotic expansions in powers of these small parameters yields an 
approximate equation for the 3 dimensional pressure known as the general 
quasi-geostrophic equation with appropriate boundary conditions. Further 
simplifying assumptions reduce the problem to the study of a 2 dimensional 
equation which describes the evolution of the temperature field on a surface 
that bounds the fluid. In the geophysical fluids literature this equation is 
known as the surface quasi-geostrophic equation. A derivation of this equa- 
tion and a discussion of its physical relevance can be found, for example, 
in Pedlosky |Pej . Salmon [S], Held at al [H PG Sj. The effects of viscosity 
are incorporated via a boundary layer analysis and a mechanism known as 
Ekman layer pumping produces the dissipative term in the 2 dimensional 
quasi-geostrophic equation. 

In the mathematical literature this 2 dimensional equation is often called 
the dissipative quasi-geostrophic equations (QG equation) with the word 
surface being omitted since the equation is 2 dimensional. This equation, 
for an unknown active scalar Q(x,t) representing the temperature on the 
boundary surface, is given by 

d t e + u ■ ve + (-A)^e = /, (1.1) 
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where U(x,t) is the velocity vector and f{x) is a given external force. The 
velocity is coupled with the temperature via a stream function $?(x,t): 

6 = (-A) 1/2 ^ = A*, (1.2) 

and 

U = = (8^% -d Xl V) = (R 2 e, -Ri&), (1.3) 

where R{ is the i th Riesz transform. Our analysis of (|1.1|) - (|1.3|) considers 
x in the 2 dimensional torus [0, 2ir] = T 2 and t G [0, oo). 

Both the non-dissipative and the dissipative QG equations have received 
much attention following seminal article of Constantin et al [CMTj . They 
observed a number of similar features between the full 3 dimensional Euler 
and Navier-Stokes equations and the much simpler QG equations in terms of 
possible formation of singularities. Recent results concerning the dissipative 
QG equations include \£Q ICCWl ICVl ICWl iDDl EE HH UM EHJ IKNV] [Ml 
IW] and references therein. 

The appropriate power f3 of the Laplacian in the derivation from the 
general 3D viscous quasi geostrophic models and Ekman boundary layer 
analysis is (3 = 1/2. Dimensionally the 2D QG equation with (3 = 1/2 is the 
analogue of the 3D Navier-Stokes equation. (5 = 1/2 is called the critical 
case. The first results concerning regularity of solutions to the dissipative 
QG equation were given in the simpler (but non-physical) subcritical case 
where j3 > 1/2: see, for example, Constantin and Wu |CW| . In the critical 
case, f3 = 1/2, Constantin, Cordoba and Wu [CCW| proved existence of 
a unique global solution evolving from any initial data that are small in 
L°° . Very recently, the smallness assumption was removed independently in 
breakthrough works of Caffarelli and Vasseur [CVj and Kiselev, Nazarov and 
Volberg |KNVj . In particular, Caffarelli and Vasseur |CVj used harmonic 
extension to establish regularity of the Leray-Hopf weak solution. On the 
other hand, Kiselev et al |KNV] proved the global well posedness of the 
critical dissipative QG equations with periodic C°° data. Their argument 
is based on a certain non-local maximum principle for a suitably chosen 
modulus of continuity. 

In the present article we consider the question of nonlinear instability of a 
steady solution of the forced critical QG equations. We note that the above 
mentioned references concern the case / = 0, but in order to ensure the 
existence of a large class of steady states we must consider the nontrivially 
forced problem. In particular, we need to reprove certain results that are 
known to hold for the unforced equations but not in the forced context, 
namely the nonlocal maximum principle of Kiselev et al [KNV| . 

The main result of this paper is that linear instability implies nonlinear 
Lyapunov instability for O, and hence U, in the function space L 2 . Such 
results connecting linear and nonlinear instability have been proven under 
certain restrictions for the 2D Euler equations, see Bardos et al |BGS| . Fried- 
lander and Vishik [VFj . and Lin (Lj. There the methods utilize a bootstrap 



NONLINEAR INSTABILITY FOR THE CRITICAL QG 



3 



technique where closure relies on the special property of conservation of 
vorticity which is valid for 2D Euler but not for 3D Euler, where the equiv- 
alent instability result is still unproven. This property cannot be utilized 
for the QG equation because the relation between the temperature and the 
stream function is not equivalent to the relation between the vorticity and 
the stream function in the 2D Euler equations. In fact this is one reason 
why it is conjectured that the QG equations might mimic possible singularity 
development in the 3D fluid equations. 

The result that linear instability implies nonlinear instability in L 2 for 
the Navier-Stokes equations in any dimension was proved in Friedlander et 
al [FPS] (see also the seminal text of Yudovich [Y]). In this case the special 
ingredient that permits the bootstrap argument to close is the smoothing 
property of the Laplacian with respect to the nonlinear term. The arguments 
in [FPSj carry over directly to the subcritical dissipative QG equation (i.e. 
(3 > 1/2) because the dissipative term again smooths the nonlinear term in 
(jl.ip - (jl.3p . However the case of the critical QG equation is more subtle 
because the critical dissipative term {(3 = 1/2) and the nonlinear term are 
now of the same order. 

Hence to prove linear instability implies nonlinear instability in L 2 for 
the critical dissipative QG equations via the bootstrap argument requires a 
different special ingredient. The one we use in this article is the existence 
of a global bound on ||V0(i)||Loo. This result for the unforced critically 
dissipative QG was proved in |KNVj and a recent preprint of Kiselev and 
Nazarov [KN] shows that the result also holds for the equation augmented 
by a dispersion term. The existence of this global bound for the forced 
equations is proven in Section [5j 

We note that the fairly general abstract theorem of Friedlander et al 
|FSV] may be applied to the critical QG equations - since the spectrum of 
the linearized operator is discrete (see Section [3]) and so the spectral gap 
condition is satisfied - and shows that linear instability implies nonlinear 
instability in H s , with s > 2. The novel result of this present paper is to 
prove instability in the "physically natural" energy space L 2 . 

Organization of the paper. In section [2] we formulate the stability prob- 
lem in terms of the temperature Q(x,t) perturbed about a steady state 
9q{x) 6 C°°. Also in the same section we define nonlinear stability /instability 
and we state the main instability result, Theorem 12.11 In section [3] we 
study the linear operator L for the dissipative QG equations in perturba- 
tion form. This operator is elliptic of order 1, with compact resolvent, and 
hence its spectrum is purely discrete for x 6 T 2 . We prove certain prop- 
erties of L that we will use in the bootstrap argument. Then in section [4] 
we use this argument to prove Theorem 12.11 In section [5] we prove, in the 
spirit of [KNV] , that the forced equation has a global C°° solution and that 
sup t>0 ||V0(t)||i°o < oo. This result in used in the bootstrap argument that 
proves the main theorem. 
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2. Notation and formulation of the result 

Let #o be the temperature of a smooth steady 2D flow with velocity qo, 
and smooth force /, that is we have 

q • W + A0 O = / (2.1) 

qo = (R20 ,-Ri6 ). (2.2) 

Here we consider 0o>9o>/ G C°°(T 2 ). We linearize (jl.ip about a the steady 
state (#o, Qo) by writing <d(x, t) = 9q(x)+9(x, t) and U(x, t) = qo(x) + q(x, t). 
In such a way we obtain an equation that governs the perturbation 9: 

d t 6 = L0 + N(6), (2.3) 

where the linear operator L is defined by 

L9 = -q -Ve-q- W - A0, (2.4) 

the velocity is coupled with the temperature via 

q=(R 2 9,-R 1 9) (2.5) 

and 

N{9) = -q-V9. (2.6) 

For simplicity of the presentation we let 9q, f, 9 have zero mean on the torus, 
and in the following we shall denote H s = {v € H S (T 2 ) : Jj 2 vdx = 0}, for 
all s > 0. We define a suitable version of stability (the same definition was 
used, e.g. in [FPS] . [VF]). 

Definition. Let (X,Z) be a pair of Banach spaces. A solution 9q of ( 12. 1D - 
(|2.2p is called (X, Z) nonlinearly stable if for any p > 0, there exists p > 
so that if 9(0) G X and \\9(0)\\z < P, then we have 

(i) there exists a global in time solution to (|2.3p such that 9(t) G 
C([0,oo);A); 

(ii) \\9(t)\\ z <p for a.e.te [0,oo). 

An equilibrium 9q that is not stable (in the above sense) is called Lyapunov 
unstable. 

The Banach space X is the space where a local existence theorem for the 
nonlinear equations is available, while Z is the space where the spectrum of 
the linear operator is analyzed, and where the instability is measured. In 
the case of the critical dissipative QG we let X be the critical Sobolev space 
H 1 (cf. [CCl [CWJ [DDI [UllJSllM]), while the growth of the perturbation is 
considered in the energy space Z = I? . Now we are ready to formulate the 
main result of the present paper. 
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Theorem 2.1. Suppose that 6q is a smooth mean-free steady state solution 
of the critical dissipative QG, i.e., it solves (]2.1j) -( i2T2l) . If the associated 
linear operator L, as defined in (|2.4j) . has spectrum in the unstable region, 
then the steady state is (H 1 ,!, 2 ) Lyapunov nonlinearly unstable. 

3. Linearized dissipative QG 

The linear operator L denned in (12. 4ft via 

L9 = -q -Ve - q- W - A# 

is a pseudo-differential operator with principal symbol 

a(x, k) = — \k\ + iqo(x) • k, 

which does not vanish on T 2 x 7? \ {0}. Therefore L is elliptic of order 1. 
Since qo, V#o €E C°°, for large enough a > 0, we have that (L — al)~ l is a 
bounded operator from L 2 into H . Moreover, the domain of L 

V{L) = {v e H l {T 2 ), f vdx = 0} C L 2 (T 2 ) (3.1) 

is compactly embedded in L 2 by Rellich's theorem, so that resolvent (L — 
a/)" 1 is a compact operator. Thus L has discrete spectrum. 

Let fi be the eigenvalue of L with maximal positive real part over L 2 . 
Let A = Re \i and </> £ L 2 be the corresponding eigenfunctionj. For a fixed 
< 5 < C\, where C\ > is a constant depending on A to be determined 
later, we denote by L$ 

L s = L-{\ + 6)1. (3.2) 

The shift ensures that L$ generates a bounded Co-semigroup over L 2 and 
that the resolvent set of L$ contains the right half plane. The following 
lemma shows that L$ generates an analytic semigroup over L 2 . 

Lemma 3.1. Over L 2 the operator L$ generates an analytic semigroup. 

The proof of the lemma modifies the proof of [Pi Theorem 7.2.7], which 
shows the analyticity of a strongly elliptic operator of order 2m over L 2 , to 
the case of the linearized QG operator, which is elliptic of order 1. 

Proof. Define the operator G via 

Gv = Av + q ■ Vv + R(v) ■ V6» + 2f3v = -Lv + 2(3v (3.3) 

where we have denoted R{v) = (R2V, —Riv) and = ||V6>o||l°°- Since qo is 
divergence-free we have that G satisfies Garding's inequality 

Re (Gv,v) > \\A 1/2 v\\ 2 L2 +(3\\v\\ 2 L2 . (3.4) 



The steady flow go = (sinma;2,0) gives an example for which the operator L has 
unstable eigenvalues over L 2 . This follows from an extension of the analysis in Friedlander 
and Shvydkoy |FS| to the dissipative equations (see also Meshalkin and Sinai IMS] ). 
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In the above estimate we also used ||i?(u)||x2 < |M|l2. Similarly, for every 
v G £>(G), we have 

|Im (Gv,v)\ < \(Gv,v)\ < \\A 1/2 v\\ 2 L2 +3{3\\vf L2 . (3.5) 

Since v is a scalar, it follows from (|3.4p and (|3.5p that the numerical range 
S(G) (cf. [Pj pp. 12]) is contained in the set 

S 4o = {X G C : -0 O < arg A < tf }, (3.6) 

where $o — arctan(3) < ir/2. Choosing $q < $ < tt/2 and defining Y,$ = 
{z £ C : argz| > we have that there is a constant C = C(#, #o) > 
such that 

dist(z,S(G)) > C|*|, for all z G S^. (3.7) 

We now claim that all real x < are in the resolvent set p{G) of the operator 
G. Recall that G = —L + 2/3/, and moreover that the spectrum of the 
operator L is contained in the half plane {z G C : Re z < A}, where 
< A = Re /x, and /x is the eigenvalue of L with largest real part with 
associated eigenfunction </>. Since qo is divergence free we also have that 

//Mia = 0) = -\\A 1/2 0\\ 2 L2 - (R(<P) ■ Wo, 4>), (3-8) 

and by taking real parts this implies that A < HV^oIIl 00 — P; hence the 
spectrum of G is contained in the right half plane, proving the claim. 

We have hence proven that is contained in the the complement of 
S(G) and has non-empty intersection with p(G); by [El Theorem 1.3.9] we 
have that C p(G) and for every z G we have the resolvent estimate 

Therefore — G is the infinitesimal generator of an analytic semigroup (cf. [Pj 
Theorem 2.5.2]) and so Lg = —G + (2/3 — A — 8)1 generates an analytic 
semigroup on L 2 , since it is a bounded perturbation of — G. □ 

Now we state and prove the lemma that will be used in the proof of our 
main result, Theorem 12.11 

Lemma 3.2. For < 7 < 1 there exists a constant C > such that 

We^vU^L* < ^HltlA- 1 ^, (3.10) 
for all smooth functions v G L 2 , where C = C(j,5,a,6o). 

Proof. Since go is divergence free, it is convenient to use the operator A a , 
defined via 

A a v = -q -Vv - Av - av = Lgv + R(v) ■ V6> - (a - A - 5)v, (3.11) 

where a > max{A + 5, C||#o||#2+e}) e > 0, and C is a sufficiently large di- 
mensional constant. We treat Lg as a bounded perturbation of A a . The 
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operator A a is also elliptic and has discrete spectrum, so by possibly choos- 
ing a different a, we have that A~ l G C{L 2 ). 
First, we claim that 

\\A- l Av\\ L 2 <C\\v\\ L 2, (3.12) 

for all smooth v G I? with zero mean. In order prove this, denote h = 
A~ l Av , which also has zero mean, and observe that (|3.12p is equivalent to 

\\h\\ I fl<C\\hr 1 A a h\\v*. (3.13) 

The definition of h implies that 

(A^M) = -(A~ 1 (q • V/i), h) - \\h\\l 2 -a\\k- l l 2 h\\l 2 , 

and therefore 

\\hf L2 + a\\A~ l l 2 h\\ 2 L2 < IIA^^II^II^II^ + Kgo-V^A- 1 ^)!. (3.14) 

Note that (go ■ VA _1 / 2 /i, A -1 / 2 /i) = since divgo = 0. Using Plancherel's 
theorem, we write this inner product in terms of Fourier coefficients (cf. |KVj 
and references therein) 

(g • V/i, A^h) = (g ■ Vh, A^h) - (g ■ VA" 1 ^, A~ l / 2 h) 

= z(2vr) 2 Y, %rk{\l\- l/2 -\k\- l / 2 )h k \l\- l,2 h. 

j+k+l=Q 

(3.15) 

In the above summation, the Fourier frequencies j, k,l G Z 2 \ {0} because go 
and h are mean free, and hk denotes the k th Fourier coefficient of h. Since 
|/| = \ j + k\ the triangle inequality gives ||/| — \k\\ < \ and therefore 

1 1 

~ W 12 

Therefore, by (I3.15P and the Cauchy-Schwartz inequality we have that 

\{q Q -Vh,A~ l h)\<C \m^\\h k \\l\- ll2 \h\ 
j+k+l=0 

<c Y b'llfol E \h-i-i\\i\- 1,2 \k\ 

jGZ 2 \{0} /gZ2\{0,-j} 

<c\\h\\ L 4A- l i 2 h\\ L , Y b'l^l^-lbr 1 -' 

j'ez 2 \{0} 

< C\\h\\ L 2\\A^ 2 h\\ L2 \\A 2 ^9 \\ L2 . (3.16) 
We plug this estimate into (|3.14p and obtain 

±\\h\\h + (a - ClIA^oll^llA-V^l^ < llA- 1 ^^ 
Since a > C\\9o\\ 2 H2+e , the above estimate proves f|3. 13|) and A~ l A G £(L 2 ). 



< ,, jKi i/2 -ifei i/2 i < \m < 

-II |||l/2|jfe|l/2 - |j|l/2|fc|l/2(| Z |l/2 + |fc|l/2) ~ 
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Now we prove that for smooth v G L 2 we have 

\\L- & l A a v\\ L 2<C\\v\\ L 2, (3.17) 

for a sufficiently large constant C > 0. The inequality (|3.17p follows by 
writing 

L^A a v = v + L^iRty) ■ W ) -(a-5- X)LJ 1 v, (3.18) 

and noting that the operator LJ 1 is bounded on L 2 (cf. [El Lemma 2.6.3]). 
Together with the boundedness of the Riesz-transforms on L 2 , (|3.18p implies 

\\Lj l A a v\\ L2 < H i2 (l + C(||W ||L°°+a-<5-A)), (3.19) 

which proves (I3~T71) and therefore Lj 1 A a G £(L 2 ). 

In order to conclude the proof of the lemma we use the fact that Lg gen- 
erates an analytic semigroup (cf. Lemma [3TTj) and therefore (cf. (Pj Theorem 
2.6.13]) we have that 

We^vWL^L* = \\L]e L ^Lpv\\L^ < (3-20) 
Now we bound ||i ( ^ 7 v||z / 2 by interpolating (cf. [Pj Theorem 2.6.10]) as follows 
WLsMl* = \\L]^{Lfv)\\ L , < CWvWl^WLfvWl, 

<C||HI^ 7 ||(^ 1 ^)(A- 1 A)(A- 1 U )||J 2 

< C||^||^ 7 ||A^||I 2J (3-21) 

where in order to obtain (|3.2ip we used (|3.17p and ()3. 12|) . Now we conclude 
the proof of the lemma by combining (|3.20p and (|3.2ip • O 

4. Proof of Theorem 12.11 

Here we prove Theorem 12.11 In order to do this we must show that the 
trivial solution 6 = of (|2.3p is (H 1 , L 2 ) Lyapunov unstable. With this goal 
in mind, we consider a family of solutions 6 s to 

d t 9 £ = L6 £ + N(6 £ ), (4.1) 

0%=o = £0, (4.2) 

where (j> is as above an eigenfunction of L associated with the eigenvalue 
with maximal positive real part A. We will prove the following proposition 
that clearly implies the desired Lyapunov instability result. 

Proposition 4.1. There exist positive constants C and e < 1 such that for 
every e G (0,e), there exists T £ > such that \\6 £ (T £ )\\ L 2 > C. 

We remark that if 9 £ (x,t) solves (|4.ip - (|4.2p . then the function Q £ (x,t) = 
9 £ (x,t) + 6>o(x) solves the forced QG equations (|5.ip - (|5.3p . with initial data 
O £ (x,0) = 6» o (x)+e0(x) G C°°(T 2 ). Moreover, in Lemma Oof Section [5] we 
prove that the global smooth solution of the forced QG equations satisfies 
||V0 e (t)||Loo < Cg for all t > 0, where the constant Cq depends solely on the 
L°° and W 1,oc norms of the initial data and the force. For e G (0, 1], we have 
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||© £ (0)||l°o < ||0o||l°° + \\<P\\l°°, and similarly ||V9 £ (0)|| L °° < ||V0 o ||l°° + 
||V(/>||l°°, which are independent of e, and therefore there exists a fixed 
Co > such that ||VG £ (t)|| LO o < C 0) for all e £ (0, 1] and for all t > 0. We 
refer the reader to the proof of Lemma 15. II for further details. The triangle 
inequality then implies that by possibly increasing Co we have 

sup||V0 e (t)|| £O o < Co (4.3) 

for all e £ (0, 1] . We will henceforth denote 9 s simply as 6 and will use the 
analogous notation for q. All constants in the following are e-independent. 

Proof of Proposition [J7} For R> := ||^||i2 to be chosen later, let T = 
T{R, e) be the maximal time such that 

||0(t)|Ua < £Re xt , for t G [0,T]. (4.4) 

Clearly T S (0, oo] due to the strong continuity in L 2 of t i— > 9(t) and the 
chosen initial condition. 

Using Duhamel's formula we write the solution of f)4. If) — (|4.2H as 

9(t) = e Lt e(/) + B(t), (4.5) 

where 

B(t) = f e L ^N(9)(s) ds. (4.6) 
First, we shall prove that 

/ . \ 1+7/2 

\\B(t)\\ L 2 <d (eRe xt ) , (4.7) 

where 7 £ (0, 1) and C\ = C(Co, \,6, 7) > are constants. To show (|4.7|) . 
we rewrite the operator B and then use Lemma 13.21 as follows: 

\\B{t)\\ L 2 = || I e {x+5) ^ s) e Ls{t - s) N{9{s)) ds\\ L 2 
Jo 

< f e ix+5){t - s) \\e Ls{t ~ s) N(9(s))\\ L 2^ L 2 ds 
Jo 

(4.8) 

where 7 £ (0, 1) is arbitrary, and C > 0. In order to bound the factor 
||A~ 1 A r (6'(s))|| i 2 we recall the explicit representation (cf. |CC] ) of the non- 
linear term 

A-\R{0) ■ V9) = C n (R 1 (9R 2 (9)) - R 2 (9R l {9))) , (4.9) 

for some dimensional constant C n , and the fact that the Riesz transforms 
are bounded on L 2 and L 4 , to obtain that 

IIA-^OO)!^ < C\\0IU0\\v < C\\9\\ 2 Li , (4.10) 
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By interpolating, we have 

11*11* < c\\otg\\e\?g . (4.ii) 

On the other hand by the Gagliardo-Nirenberg inequality and the Holder 
inequality we have that (cf. [N] ) 

11*11* = ll* 8/3 llK 8 <c||v(^/ 3 )||^ 5 

< C\\9 5 / 3 Vef/ b 8 /5 < C\\6\\% 8 \\V0f/*. (4.12) 
By combining (|4.10|) with (|4.1ip and (|4.12jl we obtain 

WA- 1 N(6(s))\\l 2 <C\\0\\% /2 \\V0\\]g (4.13) 

On the other hand, by Holder's inequality, and the boundedness of the Riesz 
transforms on L 2 , we have 

ll^)||^<||0||^||V0||^. (4.14) 

Recall that by (|4.3p we have ||V0(£)||loo < Co, for all t > 0. Using assump- 
tion (|4.4p and the fact that < 5 < C\ = A7/2, we substitute the bounds 
(I4TT31) and (ETT41) into (gSJ), to conclude 

/ s 1+7/2 

WB^W^KCAeRe^) , (4.15) 



for some positive constant C\ = C(Cq, A, 6, q, 7), proving (|4.7p . The Duhamel 
formula (|4.5p and the bound (|4.7p imply 

||0(t)|| L 2 < C^e A * + Ci (eiie^) 1+7/2 . (4.16) 



Observing that R was chosen such that R > C^, it follows that we have the 
following estimate on the maximal time T: 

«"*G^r (4 - i7> 

which clearly holds if T = 00. On the other hand, if T is finite, (|4.17p is 
obtained by combining the continuity of t 1— > \\0(t)\\ L 2, (|4.4p and (|4.16p to 
obtain 

eRe XT < C^ee XT + C 1 R 1 -*' 2 ee xr (ee AT ) 7/2 , 

which, in turn, implies (|4.17p . Therefore we have T >T e , where we defined 

1 Co 

T £ = -ln^. (4.18) 

A £ 

To conclude the proof we must find a lower bound on ||0(T £ )||^2. We use 
Duhamel's formula (|4.5p . the triangle inequality, and (|4.15p to obtain 

/ \ 1+7/2 

\\9{T e )\\ L 2 > Ctpee^ - d (eRe XT A . (4.19) 
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Using (|4.18p . with C2 given by (|4.17|) . the lower bound (14.191) implies 

R — Crk . 



= C 2 {2C 4> -R) :=C>0, 

by choosing < R < 2C^. This concludes the proof of the proposition 
which, in turn, implies Theorem l2.11 □ 



5. Global well-posedness for the forced QG equation 

In this section, by modifying the argument of Kiselev et al |KNV] . we 
prove that the forced QG equation has a unique global smooth solution. 
More precisely, we prove the following: 

Lemma 5.1. Assume that Go,/ G C°° are T 2 -periodic functions with zero 
mean. Then there exists a unique global in time smooth solution of 

d t @ + U • VG + AG = /, (5.1) 
U = R(@) = (i?2©, —Ri&), (5.2) 
6(0) = G . (5.3) 

Moreover for all t > we have 

||V9(i)|| L oo < C , (5.4) 

where Cq = CoCIIQolli 00 ||V0o||z°°, ||/||l°°, ||V/||z,«>) is a positive constant. 

The proof of the lemma is in the spirit of [IK NV] . but we additionally 
need to treat the force term, which a-priori could cause growth of the solu- 
tion. Since @(t) is mean free, it can be shown a-priori that ||G(i)||i,p, with 
2 < p < 00, remains bounded for all time. However the same methods do 
not work for the subcritical quantity ||V0(t)||z,°°, and therefore we need to 
prove the nonlocal maximum principle of jKNVj for the forced QG equation 
(|5.ip - (l5.3p . This is achieved by suitably choosing a scaling parameter B and 
making use of the fact that due to periodicity we do not need to consider 
arbitrarily large length scales. We note that the scaling parameter B is used 
only in the modulus of continuity, whereas the solutions to (15.1j) - (15.3j) are 
not space-time rescaled. 

Proof of Lemma \5.1\ We recall that a continuous, increasing, unbounded, 
concave function to : [0, 00) — > [0, 00), with w(0) = is a modulus of conti- 
nuity for a function / if 

\f(x)-f(y)\<u(\x-y\), (5.5) 

for all x,y G M. 2 . The modulus is strict if the strict inequality holds in 
(|5.5j) . We consider a modulus of continuity that also satisfies u/(0) < 00, 
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and lini£_ > o+ = — oo, namely, as in [KNV) we let 

W(0 = V-5 3 / 2 '+ 7 log(l + |log(^)), ' (5 - 6) 

where 5 > 7 > are sufficiently small fixed constants. 

Since ©o £ C°°, there exists a sufficiently large B > such that 0o has 
strict modulus of continuity = u(B£). The scaling parameter £> may 

be chosen as 

B = C||Ve ||i- exp(exp(C||e |U-)), (5.7) 

where C is a sufficiently large positive constant. Moreover, since u is un- 
bounded, by possibly increasing B we may ensure that 

AB 2 > ||V/|| L ~, (5.8) 
where the fixed dimensional constant A is as in |KNVl Lemma], and also 

^>4H|/|k-, (5-9) 

where d = diam(T 2 ) = 2tt\/2 will be fixed throughout this section. We fix 
a B that satisfies (|5.7p - (|5.9p and recall that the modulus of continuity is 
given by 

m inn \Bt-{B(f 2 , < ( < |, 

^(0=-(B0 = |,_ i3 / 2+7log(1+ . log(f)) , {>i ("0) 

Denote uj b '(0 = Bu'{B£) and w B "(() = 5V(f). We claim that w B (f) is 
preserved by the evolution f)2.3f) . so that is a global solution. We extend 
0, U, Go, / to T 2 -periodic functions on M 2 . 

The first step of the proof is to show that if 0(t) has strict modulus of 
continuity lob for t £ [0, T], then there exists r > such that 0(t) has strict 
modulus of continuity wg on i £ [0, T + r). Since ||V0(i)||z,°° < u;b'(0), we 
have that 0(i) £ C°° for all i £ [0,T], and by the local regularity theorem 
(cf. [CCl El El El) for some time r > beyond T. We must show that by 
possibly shrinking r, we have that 

\e(x,t)-Q(y,t)\<u B (\x-y\) (5.11) 

for all t £ (T, T + r) and x / y £ M 2 . 

Define the compact set K = [-27r,2vr] 2 x [-27r,27r] 2 C M 4 . Since is 
T 2 -space periodic, we have that for any (x,y) £ M 4 , with x ^ y, there exist 
(x',y') £ -KT, with x' / y', such that |x' — y'\ < \x — y\, Q(x,t) = Q(x',t), 
and Q(y,t) = Q(y',t). Because ujb is increasing, if (|5.1ip holds for all 
(x',y') £ K with x' 7^ y', then we have that for all x 7^ y £ M 2 

|0(x,t) - 0(y,i)| = |0(a/,t) - 0(y',t)\ < u B (\x' - y'\) < u B {\x - y\). 

Therefore it is sufficient to prove that there exists r > such that (|5,lip 
holds for x 7^ y, with (x, y) £ K. By assumption, there exists e > such 
that ||V0(r)||i<x> < ub'(0) — 2e, and by continuity, for small enough t we 
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have that ||V9(t)||i,°o < u B '(0) - e for all t G [T,T + r). Therefore for 
(x,y) G if, with < |x — y| = £ < p, where p < mm(5/B, e 2 /B 3 ), we have 

\@(x,t)-e(y,t)\ <e||V0(i)|U<» <£(£-e)<B£-(i?£) 3/2 = ^(£), 

for all i £ [T, T + r). On the other hand, due to the continuity in time of 
\@(x, t) — @(y, t)\, the compactness of the set {(x, y) G K : \x — y\ > p}, and 
and the fact that f)5 . 1 1 j) holds at t = T, we have that there is a sufficiently 
small r > such that (|5,lip holds for all (x, y) G K, x ^ y, and t G [T, T+r). 

The second part is to rule out the case in which there exists T > and 
x / y G K 2 such that Q(x,T)-@(y,T) =u B (\x-y\) (cf. [KNVj). Note that 
by the periodicity of G, for such x ^ y G M 2 fixed, there exist a/, y' G T 2 
such that 

u B (\x - y\) = @(x,T) - @(y,T) 

= 9(x',T) - 6(y',T) < u B (\x' - y'\) < u B (d), 

and since lo b is increasing, we must have < £ = |a; — y\ < d = diam(T 2 ). 
We conclude by showing that 

i(@(x,t)-Q(y,t))\ t=T <0, (5.12) 

contradicting the fact that the strict modulus of continuity is lost at t = T. 
In the following we suppress the time dependence of and U, since we work 
at t = T fixed. 

Since O has modulus of continuity w.b(£), we know (cf. [KNVl Lemma]) 
that U has modulus of continuity where we defined 

for some positive constant A. Then as in [KNVl Section 4] we have that 

u B (£ + h\U(x)-U(y)\)-u B (0 



\(U-Ve)(x)-(U-VQ)(y)\ < 



lim 



h^0+ h 
< \U(x) - U(y)\cu B '(0 < fi B (£) WB '(£). (5-13) 

The dissipative terms are estimated as in [KNVl Section 5], namely by the 
negative quantity 

M Bi0 = I -ii^l^-ill-^ 

+ °i + -1^-0-^(0 

Lastly, the force term is estimated using the mean value theorem 

\f(x)-f(y)\<F B (0 = {l^J h ^ °*^i> (5.15) 
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Thus, in order to conclude the proof of (I5.12|) . we must show that for all 
< £ < d, we have 

«b(Owb'(0 + Fb(0 + M B (0 < 0. (5.16) 

First we treat the case < £ < 5/B. By keeping track of B, and using 
condition (15. 8h . similar arguments as in [KNV1 Section 7] show that 

n B (£W(0 + Fb(C) + M B (0 < AB 2 £(3 + log 4~ c ) + CIIV/Hloo + 

<^(^(4 + logA)_|. (i?er i/ 2 

Since we have < B£ < 5, the above quantity is strictly negative if 5 is 
sufficiently small. Note that 5 does not depend on B. 

For the case 5/ B < £ < d, we follow the estimates in [KNV1 Section 8] to 
conclude that if 7 and 5 are sufficiently small, independent of B, then 

+ F B (0 + M B (£) < + 2||/|| £ =c - 

But B was chosen so that (|5.9|) is satisfied, i.e. 2||/||^oo < UB{d)/2ird. 
Because on [5/B, 00) the function cub(£)/£ is decreasing, for any £ G [<V-^> ^] 
we have that 2||/||l°° < to 3(d) /2nd < oo B {£/) /2ir£ > . Thus 

£ + 2||/||^ - < ^ 7 + - - - j ^ < 0, 

if 7 is sufficiently small, independent of B. Therefore (|5. 16[) holds for all < 
£ < d, and so (|5.12p is proven. Therefore the solution Q(t) exists for all time 
and has strict modulus of continuity ub, which implies that ||V0(i)||z,°° < 
w _b'(0) = B for all t > 0, concluding the proof of the lemma. □ 

Remark 5.2. We note that it is also possible to adapt the De Giorgi- 
type techniques used by Caffarelli and Vasseur |CVj to treat the forced QG 
equation. First one proves boundedness of the solution in 1? using energy 
estimates, and then the similarly to |CVl Section 2] one obtains boundedness 
(not decay) for all time of @(t) in L°° and of U(t) in BM O. The second step 
is to show that the solution is actually Holder and that it remains bounded in 
this space for all t > 0, i.e. adding a smooth force does not create additional 
difficulties. Since this is already subcritical regularity, in the third step it is 
standard to bootstrap to higher regularity and prove that the VF 1,00 norm 
of @(t) is bounded in time. 
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